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Abstract 

The chiral Potts model continues to pose particular challenges in 
statistical mechanics: it is "exactly solvable" in the sense that it sat- 
isfies the Yang-Baxter relation, but actually obtaining the solution is 
not easy. Its free energy was calculated in 1988 and the order pa- 
rameter was conjectured in full generality a year later. However, a 
derivation of that conjecture had to wait until 2005. Here we discuss 
that derivation. 



1 Introduction 

In 1970 I was in England, where my wife and I stayed for five months with 
my parents in Essex. It was largely holiday, as we were on our way back 
to Australia after two years in Boston, where I had been introduced to the 
six-vertex models and the Bethe ansatz by Elliott Lieb. 

However, I did visit Cyril Domb's group at King's College, London, and it 
was there that I first interacted with Tony Guttmann, who was also visiting 
the department: he was an invaluable aid to navigating the labyrinthine cor- 
ridors and staircases that linked the department's quarters in Surrey Street 
with the main part of the College. 

Tony's natural enthusiasm for statistical mechanics must have been in- 
fectious, for it was at this time that I realised that the transfer matrices of 
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the six-vertex model commuted - a vital first step in the subsequent solution 
of the eight- vertex model. 

This led to the solution of a number of other two-dimensional lattice 
models. One that has proved particularly challenging is the chiral Potts 
model. Here I wish to discuss some of the insights that led to the recent 
derivation of its order parameters. 

The chiral Potts model is a two-dimensional classical lattice model in 
statistical mechanics, where spins live on sites of a lattice and each spin takes 
N values 0, 1, . . . , N — 1, and adjacent spins interact with Boltzmann weight 
functions W, W. We consider only the case when the model is "solvable" , by 
which we mean that W, W satisfy the star-triangle ( "Yang-Baxter" ) relations 
P]. The free energy of the infinite lattice was first obtained in 1988 by using 
the invariance properties of the free energy and its derivatives. |2] Then in 
1990 the functional transfer matrix relations of Bazhanov and Stroganov [3] 
were used to calculate the free energy more explicitly as a double integral. |U 
IHJ Ej The model has a critical temperature, below which the system exhibits 
ferromagnetic order. 

The next step was to calculate the order parameters Mi,... , M.n-i (de- 
fined below). These depend on a constant k which decreases from one to zero 
as the temperature increases from zero to criticality. In 1989 Albertini et al 
jZj made the elegant conjecture, based on the available series expansions, 
that 

Mr = k r{N ' r)/N2 , < r < N . (1) 

It might have been expected that a proof of such a simple formula would not 
have been long in coming, but in fact it proved to be a remarkably difficult 
problem. Order parameters (spontaneous magnetizations) are notoriously 
more difficult to calculate than free energies. For the Ising model (to which 
the chiral Potts model reduces when N = 2), the free energy was calculated 
by Onsager in 1944 8 , but it was five years later when at a conference in 
Florence he announced his result for the spontaneous magnetization, and not 
till 1952 that the first published proof was given by Yang P ITU]. 

Similarly, the free energy of the eight-vertex model was calculated in 
1971. [UJ The spontaneous magnetization and polarization were conjectured 
in 1973 and 1974, respectively [T^ HH] . but it was not till 1982 that a proof of 
the first of these conjectures were published [TJ]. A proof of the second had 
to wait until 1993[J3]! 

By then three separate methods had been used. The Onsager- Yang cal- 
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dilation was based on the particular free-fermion/spinor/pfaffian/Clifford 
algebra structure of the Ising model[in|. As far as the auther is aware, this 
has never been extended to the other models: it would be very significant if 
it could be. 

The eight-vertex and subsequent hard-hexagon calculation was made us- 
ing the corner transfer matrix method, which had been discovered in 1 9 76 [T7] . 
This worked readily for the magnetization (a single-site correlation), but not 
for the polarization (a single-edge correlation). This problem was remedied 
by the "broken rapidity line" technique discovered by Jimbo et al [To] . 

For all the two-dimensional solvable models, the Boltzmann weight func- 
tions W, W depend on parameters p and q. These parameters are known as 
rapidities and are associated with lines (the dotted lines of Figure [T]) that 
run through the midpoints of the edges of the lattice. In general these are 
complex numbers, or sets of related complex numbers. In all of the models 
we have mentioned, with the notable exception of the N > 2 chiral Potts 
model, these parameters can be chosen so that W, W depend only on the 
rapidity difference (spectral parameter) p — q. 

This property seems to be an essential element in the corner transfer 
matrix method: the star-triangle relation ensures that the corner transfer 
matrices factor, but the difference property is then needed to show that the 
factors commute with one another and are exponentials in the rapidities. The 
difference property is not possessed by the N > 3 chiral Potts model and one 
is unable to proceed. At first the author thought this would prove to be 
merely a technical complication and embarked on a low-temperature numer- 
ical calculation JH] in the hope this would reveal the kind of simplifications 
that happen with the other models. This hope was not realised. 

I then looked at the technique of Jimbo et al and in 1998 applied it to the 
chiral Potts model. One could write down functional relations satisfied by the 
generalized order parameter ratio function G pq (r), and for N = 2 these were 
sufficient (together with an assumed but very plausible analyticity property) 
to solve the problem. However, for N > 2 there was still a difficulty. Then 
p, q are points on an algebraic curve of genus > 1 and there is no obvious 
uniformizing substitution. The functional relations themselves do not define 
G pq (r): one needs some additional analyticity information, and that seems 
hard to come by. 

The calculation of the free energy of the chiral Potts model jSJ |HJ E] 
proceeds in two stages. First one considers a related "r 2 (t g )" model. |20. 
This is intimately connected with the superintegrable case of the chiral Potts 
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model. [21] It is much simpler than the chiral Potts model in that its Boltz- 
mann weights depend on the horizontal rapidity q only via a single parameter 
t q , and are linear in t q . Its row-to-row transfer matrix is the product of two 
chiral Potts transfer matrices, one with horizontal rapidity q, the other with 
a related rapidity r = VRq defined by eqn. (jHJ) of section 2. 

For a finite lattice, the partition function Z of the T2(t q ) model is therefore 
a polynomial in t p . The free energy is the logarithm of Z l / M , where M is the 
number of sites of the lattice, evaluated in the thermodynamic limit when 
the lattice becomes infinitely big. This limiting function of course may have 
singularities in the complex t q plane. A priori, one might expect it to have N 
branch cuts, each running though one of the TV roots of unity. However, one 
can argue that in fact it only has one such cut. As a result the free energy 
(i.e. the maximum eigenvalue of the transfer matrix) can be calculated by a 
Wiener-Hopf factorization. 

The second stage is to factor this free energy to obtain that of the chiral 
Potts model. 

It was not until 2004 that I realised that : 

(1) If one takes p, q to be related by eqn. (|T3|) below, then G pq {r) can 
be expressed in terms of partition functions that involve p, q only via the 
Boltzmann weights of the T2(t p i) model, with p' = i? _1 p. 

(2) It is not necessary to obtain G pq (r) for arbitrary p and q. To verify 
the conjecture (JTJ) it is sufficient to obtain it under the restriction (|14p. 

I indicate the working in the following sections: a fuller account is given 
in Ref. [22J. The calculation of G pq (r) for general p, q remains an unsolved 
problem: still interesting, but not necessary for the derivation of the order 
parameters M. r . 

2 Chiral Potts model 

We use the notation of PQI3JI2B]- Let k, k' be two real variables in the range 
(0, 1), satisfying 

k 2 + k' 2 = 1 . (2) 
Consider four parameters x p , y p , fi p , t p satisfying the relations 

kx N p = 1 - k'/n» , ktf = 1 - h>y» , t p = x p y p . (3) 

Let p denote the set {x p , y p , fi p , t p }. Similarly, let q denote the set {x q , y q , fi q , t q }. 
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We call p and q "rapidity" variables. Each has one free parameter and is a 
point on an algebraic curve. 

Define Boltzmann weight functions W pq (n), W pq {n) by 

w pq (n) = ^ p /^ q r ft , (*o 

W pq (n) = (w q Y ft ^ ~ ^ , (4b) 

where 

They satisfy the periodicity conditions 

W pq (n + N) = W pq (n) , W pq (n + N) = W pq (n) . 

Now consider the square lattice £, drawn diagonally as in Figure [TJ with 
a total of M sites. On each site i place a spin ai, which can take any one of 
the N values 0, 1, . . . , iV — 1. 

The solid lines in Figure [T] are the edges of C. Through each such edge 
there pass two dotted or broken lines - a vertical line denoted v and a hor- 
izontal line denoted h (or p or q). These v,h,p,q are rapidity variables, as 
defined above. We refer to each dotted line as a "rapidity line" . 

With each SW - NE edge (with i below j) associate an edge weight 
W v h(<Ji — Vj)- Similarly, with each SW - NE edge (j, k) (j below k), associate 
an edge weight W v h{oj — a^). (Replace h by p or q for the broken left and 
right half-lines.) Then the partition function is 

2 = EII W vh (ai - a 3 ) \{W vh {a j - a k ) , (5) 

a 

the products being over all edges of each type, and the sum over all N M 
values of the M spins. We expect the partition function per site 

K = Z X ' M 

to tend to a unique limit as the lattice becomes large in both directions. 

Let a be a spin on a site near the centre of the lattice, as in the figure, 
and r be any integer. Then the thermodynamic average of uo ra is 

F pq {r) = (u ra ) = Z-^^U^hi^-^UWvh^j-^k) ■ (6) 
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Figure 1 : The square lattice (solid lines, drawn diagonally) , and the associated rapidity 
lines (broken or dotted). 

We expect this to also tend to a limit as the lattice becomes large. 

We could allow each vertical (horizontal) rapidity line a to have a different 
rapidity v a {hp). If an edge of £ lies on lines with rapidities v a , hp, then the 
Boltzmann weight function of that edge is to be taken as W v h(n) or W v h(n), 
with v = v a and h = hp. 

The weight functions W pq (n), W pg (n) satisfy the star- triangle relation. [I] 
For this reason we are free to move the rapidity lines around in the plane, 
in particular to interchange two vertical or two horizontal rapidity lines. 
So long as no rapidity line crosses the site with spin a while making such 
rearrangements, the average (u ra ) is unchangedby the rearrangement. 1 

1 Subject to boundary conditions: here we are primarily interested in the infinite lattice, 
where we expect the boundary conditions to have no effect on the rearrangements we 
consider. 
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All of the v, h rapidity lines shown in Figure ^ axe "full" , in the sense that 
they extend without break from one boundary to another. We can move any 
such line away from the central site to infinity, where we do not expect it to 
contribute to (uj ra ). Hence in the infinite lattice limit F pq (r) = (uj ra ) must 
be independent of all the full-line v and h rapidities. 

The horizontal rapidity line immediately below a has different rapidity 
variables p, q on the left and the right of the break below a. This means that 
we cannot use the star-triangle relation to move it away from a. 

It follows that F pq (r) will in general depend on p and q, as well as on 
the " universal" constants k or k' . We are particularly interested in the case 
when q = p. Then the p, q line is not broken, it can be removed to infinity, 
so 

M r = F pp {r) = (uj ra ) = independent of p . (7) 

These are the desired order parameters of the chiral Potts model, studied 
by Albertini et al. By using this "broken rapidity line" approach, I was finally 
ably to verify their conjecture (0) in 2005(23 122] • Here I shall present some 
of the observations that enabled me to do this. 



Automorphisms 

There are various automorphisms that change x p , y p fi p , t p while leaving the 

relations (0 ) still satisfied. Four that we shall use are R, S, M, V, defined 
by: 

{%Rp, VRpi ^Rpi tRp} — {y P ,UX p ,l/ fi p ,U}t p } , 

{xsp,ysp,Hsp,ts P } = {l/y P A/xp,u~ 1/2 yp/(x p fi p ),l/t p } , (8) 

\%Mpi VMpi I^Mpt tMp} = \ x piyp-i^^pitp\ i 

\%vpi Vvpi l^vp: tvp} — \ x pi vyp, fl p , Ut p } . 



The central sheet V and its neighbours. 

We shall find it natural, at least for the special case discussed below, to 
regard t p as the independent variable, and x p , y p , \i p to be defined it terms of 
it by (JsJJ). They are not single- valued functions of t p : to make them single- 
valued we must introduce N branch cuts B , B\, . . . , Sjv-i in the complex 
tp-plane as indicated in Figure (J2J. They are about the points 1, u, . . . , uj 1 ^^ 1 , 
respectively, 
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Figure 2: The cut t p -plane for N = 3. 

Since the Boltzmann weights are rational functions of x p ,y p , we expect 
G pq (r), considered as a function of t p or t q , to also have these N branch cuts. 

Given t p in the cut plane of Figure 121 choose fi p to be outside the unit 
circle. Then x p must lie in one of N disjoint regions centred on the points 
1,uj, ... , u 1 ^' 1 . Choose it to be in the region centred on 1. We then say that 
p lies in the domain T>. When this is so (and t p is not close to a branch cut), 
then in the limit k' — > 0, n p = 0{l/k') and x p — > 1. 

The domain V has N neighbours Vq, . ..,Pjv_i , corresponding to t p 
crossing the N branch cuts Bo, . . . , B^-i, respectively. The automorphism 
that takes V to Vi, while leaving t p unchanged, is 

Ai = V l - 1 RV N ~ i . (9) 

The mappings A4 are involutions: A\ — 1. 
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3 Functional relations 



We define the ratio function 



G Pq (r) 



F m (r)/F pq (r-1) . 



(10) 



The functions F pq (r), 



G Pq ( r ) 



satisfy two reflection symmetry relations. 
Also, although we cannot move the break in the (p, q) rapidity line away 
from the spin a, we can rotate its parts about a and then cross them over. 
As we show in and , this leads to functional relations for G pq {r): 



G 



Rp,RqV 

G p>g (r 
G pq {r 

Gpqir 



GMp,q(r 
N 

II G P q(r 

r=l 



l/G„(N-r + l) , 
l/G RSq , RS p(N-r + l) 
GR q: R-ip(r) , 

•KqfJ'q ^ ^pf^p 



(11) 



Vp^q 



G 



p,M~ 



Vqf^p 

= G 



Gr-^^r p { 
r + 1) , 



Also, from ((7|), 

M r = G pp {l)---G pp {r) . (12) 

For the case when N = 2 we regain the Ising model. As is shown in 
there is then a uniformizing substitution such that x p ,y p , fi p ,t p are all 
single-valued meromorphic functions of a variable u p , and W pq (n), W pq (n) 



u p . In fact 



and hence G pq (r) depend on u p , u q only via their difference w, 
all quantities are Jacobi elliptic functions of u p , u q with modulus k. One can 
argue (based on low-temperature series expansions) that G pq {r) is analytic 
and non-zero in a particular vertical strip in the complex u q — u p plane. The 
relations (|llj) then define G pq (r). They can be solved by Fourier transforms 
and one readily obtains the famous Onsager result 



k' 2 ) 1/s 



(13) 



For iV > the problem is much more difficult. There then appears to be 
no uniformizing substitution and G pq (r) lives on a many-sheeted Riemann 
surface obtainable from T> by repeated crossings of the branch cuts. One 
can argue from the physical cases (when the Boltzmann weights are real and 



9 



positive) that G pq (r) should be analytic and non-zero when p,q both lie in 
V, but the relations only relate these sheets to a small sub-set of all 
possible sheets. There seems to be a basic lack of information. 

4 Solvable special case: q = Vp 

The author spent much time mulling over this problem, then towards the 
end of 2004 he realised that the case 

q = Vp (14) 

may be much simpler to handle, and still be sufficient to obtain the order 
parameters M. r . 

The reason it is simpler is that one can rotate the left-half line p anti- 
clockwise below a until it lies immediately below the half-line q, as in Fig. 

5 of [22]- One has to reverse the direction of the arrow, which means the 
rapidity is not p but p' = R^p. 

The result is that p enters the sums in (jSJ), © only via the weights of the 
edges shown in Figure El The left-hand spins are the same - the spin a. The 
right-hand spins are set to the boundary value of zero. 

Further, we can sum over the spins between lines p' and q. For instance, 
summing over the spin g gives a contribution 

U(b,c,d,e) = J2 W vp'( b ~ gWv P '(c- g)W vq (g - d)W vq (g - e) . 

9 

If a, a%, . . . , (Tl are the spins on the lowest row of Figure 01 and a, a[, . . . , a' L 
are those in the upper, then the combined weight of the edges shown in 
Figure El is 

n%-i,<T i)ff >u) . (is) 

i=i 

Now q = VRp', which from (jHJ) means that 

x q = y P ' , y q = uj 2 x p i , [i q = l/7y • (16) 

This is the equation (3.13) of [3], the q, r therein being our p', q and k, £ having 
the values 0,2. From (3.17) therein, U(b,c,d,e) vanishes if < mod(6 — 
e,N) < 1 and 2 < mod(c — d, N) < N — 1. It follows that the spins in 
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Figure 3: The lattice after rotating the half- line p to a position immediately below q. 

the upper row are either equal to the corresponding spins in the lower row, 
or just one less than them. From (2.29) and (3.39) of it follows that to 
within "gauge factors" (i.e. factors that cancel out of eqn. ITHJl U (b, c, d, e) 
depends on p very simply: it is linear in t p . 

In fact, these Boltzmann weights U(b,c,d,e) are those of the T 2 (t p i) 
modeljH E| mentioned earlier. Just as this model plays a central role 
in the calculation of the chiral Potts free energy, so it naturally enters this 
calculation of the order parameters. 

In the low-temperature limit, when k' — > 0, fx g ~ 0(k'~ 1 ^ N ), x p , x q — > 
1, we can verify that the dominant contribution to the sums in (jSj), © comes 
from the case when a±, . . . , <7l, a[, . . . , a' L are all zero. Also, to within factors 
that cancel out of lfTKjl and (0), 

U(b, c, c, b) = 1 — utpi = 1 — t p . (17) 
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It follows that the RHS of (JHJ), and therefore of (jlUJl . is a ratio of two 
polynomials in t p , each of degree L, and each equal to (1 — t p ) L in the limit 
k' —>■ 0. By continuity (keeping L finite), for small values of k' their L zeros 
must be close to one. Provided this remains true (which we believe it does) 
when we take the limit L — > oo, we expect G p y p {r) to be an analytic and 
non-zero function of t p , except in some region near t p — 1. As k! becomes 
small, this region must shrink down to the point t p = 1. 

Similarly, if we rotate the half line p in Figure ^ clockwise above a, we 
can move it be immediately above q, with p replaced by Rp, as in Fig. 6 of 

• The p', q of Figur^21 herein are now replaced by q, Rp. This corresponds 
equation (3.13) of [4 with the q,r therein replaced by q,Rp. From (fTl|) 
it follows that k,£ in |3] now have the values —1,N + 1. The combined 
star weights U are now those of the Tjv(^p) model. They are polynomials in 
t p of degree N — 1, except for terms which contribute a factor x p ^ to the 
contribution of (fT5j) to G p y p (r), where 

e(r) = l-N6 r ,o , (18) 

the 5 function being interpreted modulo N, so e(0) = e(N) = 1 — N. 

When k! — > these polynomials are (1 — ut p )(l — uj 2 t p ) ■ ■ ■ (1 — u N ~ l t p ). 
In the large-L limit, with k' not too large, we therefore expect x p ^G p y p (r) 
to have singularities near t p — u, . . . , lj 1 ^' 1 , but not near t p — 1. 

If we define 

g(p;r) = G p y p {r) , (19) 

then this implies that the function x e p ^g(p; r) does not have Bq as a branch 
cut. This is in agreement with the fourth and sixth functional relations in 
(|11|) . If we set q = Vp therein we obtain 

x;^g(p;r) = y^giV^Rp-r) , (20) 

using V~ l R = R~ l V . Here we have used the fourth relation for r / 
and the sixth to then determine the behaviour for r = 0. (For r = the 
fourth relation merely gives = 0.) From (jUJ) the automorphism V~ l R is 
the automorphism A that takes p across the branch cut B , returning t p 
to its original value, while interchanging x p with y p . Thus (j2l?j) states that 
x p t ^ s> g(p]r) is the same on both sides of the cut, i.e. it does not have the 
cut B . 

These are the key analyticity properties that we need to calculate g(p; r) 
and Ai r . We do this in [221 I2H]; but this meeting is in honour of Tony 



12 



Guttmann, an expert in series expansion methods, so it seems appropriate 
to here describe the series expansion checks I made (for N = 3) when I first 
began to suspect these properties. 



5 Consequences of this analyticity 

The above observations imply that g(p;r), considered as a function of t p , 
does not have the branch cuts of Figure |2l except for the branch cut on the 
positive real axis. 

This means that g(p; r) is unchanged by taking allowing t p to cross any 
of the branch cuts Bi, . . . , -B/v-i and then returning it to its original value, 
i.e. it satisfies the N — 1 symmetry relations: 

g(p;r) = g(Aip;r) for % = 1, . . . , N - 1 , (21) 

Ai being the automorphism (JSJ). 

For N = 3, this can be checked using the series expansions obtained in 
[213 • We use the hyperelliptic parametrisation introduced in (2Z1 EH EH] • We 

define parameters x,z p ,w p related to one another and to t p by 



n 3n ' 



W 



n 7T 

n=l I 1 



(^) 2 = 27,n(j-^j • (22) 
(1 - x 2n - x zjw){\ - x 2n - l w/z){l - x 6n ~ 5 zw)(l - x &n ~ l z- l w- 1 ) 



x 2n 2 z /w)(l — X 2n w/z)(l — X %n 2 ZVj)(l — X 6n A Z 1 W 1 

(writing z p ,w p here simply as z,w), and 



(23) 



f(uz p ) f(-u/w p ) 2 f(-UW p /z p ) 

lp f(u 2 z p ) f{-u 2 /w p ) f(-uJ 2 w p /z p ) ' 1 } 

where f(z) is the function 

oo 

/(*) = \{{l-x n - l z){l-x n /z) . (25) 

n=l 

Note that x, like k', is a constant (not a rapidity variable) and is small 
at low temperatures. We develop expansions in powers of x. For p in T>, the 
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parameters Zpj ti)p ctrG of order unity, so to leading order w p = z p + 1, x p — 1, 
y p = (uj- u?z v )j (1 - u; 2 z p ). 

The automorphisms R, S, V transform z p , w p to 

zr p = xz p , z Sp = l/(xz p ) , z Vp = —l/w p 

w Rp = z v /w p , w Sp = l/(xw p ) , w Vp = z p /w p , (26) 
so from @, if p% = Aip then 

Z po 1/ (OTj,), Z pi XWpj Zp, Z p2 Zp 

w po = -l/(xz p ), w pi = Wp, w P2 = xzp/vjp . (27) 

If we write g{p;r) more explicitly as g(z p ,w p ;r), then the relations (J2"T|) be- 
come 

g(z p ,w p ;r) = g(-xw p /z p ,w p ;r) (28a) 
g(z p ,w p ;r) = g(z p ,xz p /w p ;r) . (28b) 

Using (j5J), (jHJ), we can write (fTU|) as 

G Pg (r) = E^F M (j) lj^u^F pq {j) , (29) 

3=0 I j=0 

where F pq (j) is the probability that spin a has value j. 

We use the series expansions (39) - (52) of [26J for F pq (l) / F pq (0) and 
F pq (2)/F pq (0) in terms of 

a = z q /zp , /3 = Wg/u^ . (30) 

Since g = Up, z g = —l/w p , w q = z p /w p and we find from (39) of that 
u = —uJWp/zp. (Choosing the cube root for u to ensure that F pq (i) / F pq (0) is 
real when y p = y q = which is when z p = uj 2 , w p = —uj: we then regain the 
physically interesting q = p case of eqn. [7| ) For p, q in T>, the parameters 
z p ,Wp, z q ,w q ,a, (3 are all of order unity, we can then use the expansion (48) 
of j2S] to obtain 

F pq (l) / F pq (0) = uj 2 ^x{z v ) = uj 2 i) 2 {-Wp) , 
F pq (2)/F pq (0)=ujilj 2 (z p ) = uj^i-Wp) , (31) 
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where 



^(z) = -( z + l)x + (z + lfx 2 /z - (z 3 + 6z 2 + 16z + 16 + Az' 1 + z~ 2 )x 3 
+ (z A + llz 3 + Alz 2 + 85^ + 81 + 25z _1 + lz~ 2 + ^ 3 )x 4 + 0(x 5 ) , 

and 

il> 2 (z) = zx - (2z + 1 + z~ l )x 2 - (z 2 - 8z - 2 - 3z' 1 - Z - 2 )x 3 

-(2z 3 - 5z 2 + 31 z + 6 + Uz- 1 + hz~ 2 + z~ 3 )x A + 0(x 5 ) . 

The automorphism ()28a|) interchanges "D with T>\. To leading order in x, 
the mid-point is when z p = ix 1 ^ 2 ,w p = 1. This is on the boundary of the 
domain T>, in which the series (48) of was obtained, so the series is not 
necessarily convergent at this point. Nevertheless, if we take z p = 0(x 1 ^ 2 ) 
in the above two series, we find the terms originally of order x^ become of 
order not larger than x^ +1 ^ 2 . Extrapolating, this suggests that the series 
do still converge at the midpoint, so we can use them to check whether the 
symmetry is satisfied. 

The first check occurs at order x 3//2 , where both series contain a term 

± (xz p — x 2 w p /z p ) 

(using the fact that to leading order w p = 1 at the midpoint). This is indeed 
symmetric under z p — > —xw p /z p . If we subtract this term from the series 
(using the expansion of w p in terms of z p ), we can then check the behaviour 
at order x 2 , and similarly then at order x 5//2 . All three checks are satisfied 
by both series. 

The perceptive reader will remark that (pHj) allows us to work with w p 
instead of z p . Since w p is unchanged by Ax, the symmetry appears obvious. 
Indeed it is, but only because a quite remarkable event occurred in deriving 
these series, namely the z series contains no powers of z + 1 as denominators, 
and the w series contains no powers of w — 1. If one expands w in terms of 
z (or z in terms of w), then one does find such terms. It is their absence 
from (|3ip that makes the series obviously convergent near w = 1 or z = — 1. 
I have presented the argument in terms of z p to make it clear that one does 
indeed have three non-trivial checks on the symmetry to the available order 
of the series expansion. 
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Similarly, (|28b|) interchanges D with D 2 , with mid-point z p = —l,w p = 
ix 1//2 . If one now works with w p as the variable, one can verify to the same 
three orders the symmetry w p — > xz p /w p . 

So our series provide no less than six checks on the symmetries ()28a|) . 
f!28b|) . When I first observed this, I could see the resemblance to the prop- 
erties of the free energy of the T2(t q ) model. One such property is that 
T2(t q )r2(ujt q ) ■ ■ ■ r 2 (uj N ~ 1 t q ) is a rational function of x^ , so I looked at the 
series for 

N-l 

C(p;r) = \{g{V ] p]r) 

j=0 

= g(z p ,w p ;r) g(-l/w p , z p /w p ;r) g(-w p /z p ,-l/z p ;r) . (32) 

Choosing an arbitrary value for z p and working to 30 digits of accuracy, I 
soon found that the series (known to order x 4 ) fitted with the simple formulae 

£(p;0) = l/x 2 p , C(p; 1) = k l/3 x p , £(p; 2) = k- 1/3 x p . (33) 

All this strongly suggested that I was on the right track. It did not take 
long to justify my observations for general N. For instance, if g(p; r) only has 
the branch cut B , and x~ e ^g(p; r) does not have that cut, then x~ e ^C(p; r) 
does not have the cut B . But this function is unchanged by p — > Vp, which 
rotates the t p plane through an angle 27r/N. Hence it cannot have any of the 
cuts B Q , Bi, . . . , Bjsr-\. We do not expect any other singularities (e.g. poles) 
for p in V, so the function is analytic in the entire t p plane. It is bounded 
(the Boltzmann weights W, W remain finite and non-zero as y p — * oo, the 
ratio Hp/ y p remaining finite), so from Liouville's theorem it is a constant 
(independent of p but dependent on r). 

We can relate these constants to the desired order parameters Ai r in two 
ways, and then use these relations to calculate the M. r . When y p = y q = 
and x p = k 1 ^ , our special case q = Vp intersects with physically interesting 
case q = p, so from ((7j) , 

Xp- e{r) £(p; r) = k-< r) ' N (Mr/M r -i) N . (34) 

When y p = y q = oo (Hp/y p remaining finite) and x p = fc _1//7V we find not 
q = p but q = M~ l p, which is related to q = p by the fifth of the functional 
relations (JTTJ), giving 

x-^C(p-r) = k e ^ N {M r+1 /M r ) N . (35) 
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The left-hand sides of these last two equations, being constants, are the 
same in both equations. We can therefore equate the two right-hand sides, 
for r — 1, . . . , N — 1. Using the fact that Aio = A4n = 1, we can solve for 
Aii, . . . ,A4n-i to obtain 

Mr = k r{N ~ r)/N2 for r = 0, . . . , N , (36) 

which verifies the conjecture (JTJ) of Albertini et al [7j. For N = 3 these results 
do of course agree with my original conjectures (f3~3l . 

In j22] I also show that one can calculate Gpy p (r) = g(p\ r) by a Wiener- 
Hopf factorization, giving 

g(p- r ) = k (N + l-2r)/N* S e(r) (3?) 

for r = 1, . . . , N, where 

log Sp = - W2 lo g k+— I loglAM -t p]M , (38) 

and 

A(6) = [(l-2k'cos9 + k' 2 )/k 2 f N . (39) 

(This function S p should not be confused with the automorphism S defined 
in ©. 

As is implied by the above equations, S p satisfies the product relation 

SpSvp ■ ■ ■ Svn-i p = k~ x l N x v . (40) 

Also, if one sets q = Vp in the second of the relations (fTTj) . uses the identity 
RS = MVRSV and the fifth relation, one obtains g(p; r)g(RSVp; N—r) = 1, 
from which we can deduce the symmetry 

S p S RSVp = k- 2 / N2 . (41) 

For N = 3 the automorphism p — > RSVp takes z p , w p to — w p , —z p , so 
this relation can then be written 

S(z p ,w p )S(-w p ,-z p ) = k~ 2 / g . (42) 
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6 Another interesting case: q = V 2 p 

We now have the solution for G pq (r) for q = p and for q = Vp. This suggests 
looking at one more case: q = V 2 p, where y q = w 2 y p . Similarly to section 5, 
we set g 2 (p; r) = G pq (r) and 

AT-l 

L 2 (p;r) = JJ g 2 (V 1 p;r) . 
j=o 

For N = 3 we have used the series expansions of to obtain for this 

case 

F pq {\) = u<j){w p ) , F pq {2) = tu 2 0(l/w p ) , (43) 

where 

(j>(w) = (w- l)x - (2w 2 -2w + l)x 2 /w + (2w 3 + 6w 2 - Qw + l)x 3 /w 

- (2w i + 8w 3 + 2Aw 2 - 22w + 5)x 4 /w + 0(x 5 ) . (44) 

As in the previous case, the coefficients are Laurent polynomials in w. There 
is no sign of any singularity near w p = 1, t p = uj so this suggests that G pq (r), 
considered as a function of t p , does not have the branch cut B\. 

Indeed, this is a consequence of the third functional relation (JTTJ). Setting 
q = V 2 p therein, we obtain 

92(p;r) = g 2 (A 1 p;r) , 

which tells us that g 2 {p] r) is unchanged by taking t p across the branch cut 
B\ and returning it to its original value. This means that the cut Bi is 
unnecessary. However, g 2 (p;r) does appear to have the other two cuts B 
and B 2 . 

To the available four terms in the series expansion we found 

L 2 (p;l) = x\ , 

and 

L 2 (p;0) = k- l /*x- l h{z p ,w p f , L 2 (p;2) = k l ^x- p l h(z p} w p )- z , (45) 
where 

h(z, w) — 1 + (x 2 — 6x 3 + 35x 4 )(w/z 2 + zw — z/w 2 + 3) 
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+ x\w 2 jz A + z 2 /w 4 + z 2 w 2 - 3) + 0(x 5 ) . (46) 

The result for L 2 (p; 1) looks encouraging, and indeed to the four available 
terms in the series expansion we also find 

g 2 (p-l) = k 2 ' 9 S p S Vp . (47) 

The results for L 2 (p;0) and L 2 (p;2) are not so encouraging and I have 
failed to find any obvious result for these or for g 2 (p;0), g 2 (p;2). In [22] I 
conjecture that for general N the functions G p yi p (r) have a simple form as 
a product of S functions provided i — 0, . . . , N — 1 and r = 1, . . . , N — i. For 
other values of i, r they remain a puzzle. (Except when i = 1 and r = N: 
this case can be deduced from the sixth relation of eqnlTTl) 

If (14 7|) is correct, then we have some information on the function L pq (r) 
of eqn. 56 of [23] . From this and the first equation of (fTTj). 

L pq (r) = G pq (r)G Rq , Rp (r) = G pq (r)/G qp (N - r + 1) . (48) 

Setting q = Vp and using (JT^), we obtain 

L pq (r)=g(p;r)/g 2 (Vp;N-r + l) . (49) 

Taking r = 0, it follows from (J37|) and (J47|) that 

L pq {0) = k-*/ g /(S 2 S VP S V 2 P ) = k~ l ' 9 /(x p S p ) . (50) 

The function L pq , for arbitrary p, q, was introduced in [2H] partly because its 
square is a rational function of x p , y p , fi p , x q , y q , /i g when iV = 2, so the hope 
was that it might be similarly simple for all N. We see that this cannot be 
so: S p is not such a function. 

7 Summary 

I have outlined the recent derivation of the order parameters of the solvable 
chiral Potts model, a derivation that verifies a long-standing and elegant 
conjecture. [7] As with all the calculations on solvable models satisfying the 
star-triangle relations, the trick is to generalize the model to a point where 
one has a function, here G pq (r), to calculate, rather than a constant, as 
one can obtain relations and properties that define this function. On the 
other hand, this is an example where it pays not to over-generalize: we can 
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handle the particular function G p y p (r), and this is sufficient for the purpose 
of obtaining the order parameters. The general G pq (r) continues to defy 
calculation. 

Series expansion methods can provide a valuable check on such deriva- 
tions, which are of their nature believable but hard to make fully mathe- 
matically rigorous. One usually tries to present the argument in as logical a 
manner as possible, but this is usually not the manner in which it was origi- 
nally developed. Here I have indicated the points in the calculation when I 
found the available checks both reassuring and encouraging. 
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